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ABSTRACT 

We present the results of SPH simulations in which two clouds, each having mass 
Mq = 500 Mq and radius = 2 pc, collide head-on at relative velocities of = 
2.4, 2.8, 3.2, 3.6 and 4.0kms“^. There is a clear trend with increasing At;^. At low 
A-Pq, star formation starts later, and the shock-compressed layer breaks up into an 
array of predominantly radial filaments; stars condense out of these filaments and fall, 
together with residual gas, towards the centre of the layer, to form a single large-A 
cluster, which then evolves by competitive accretion, producing one or two very mas¬ 
sive protostars and a diaspora of ejected (mainly low-mass) protostars; the pattern of 
filaments is reminiscent of the hub and spokes systems identified recently by observers. 
At high Apq , star formation occurs sooner and the shock-compressed layer breaks up 
into a network of filaments; the pattern of filaments here is more like a spider’s web, 
with several small-A clusters forming independently of one another, in cores at the 
intersections of filaments, and since each core only spawns a small number of proto¬ 
stars, there are fewer ejections of protostars. As the relative velocity is increased, the 
mean protostellar mass increases, but the maximum protostellar mass and the width 
of the mass function both decrease. We use a Minimal Spanning Tree to analyse the 
spatial distributions of protostars formed at different relative velocities. 
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1 INTRODUCTION 


It is commonly assumed that collisions between interstellar 
clouds trigger star formation, by producing dense, gravita¬ 
tionally unstable layers of gas. However, it is difficult to 
evaluate the contribution that collisions make to the overall 
galactic star formation rate, as we do not know how fre¬ 
quently such collisions occur, nor how star formation pro¬ 
ceeds when they do (e.g Pongracic et al.||1992[ Chapman 


et al. II 1992| [Whitworth et al.||1994| [Heitsch et al. |2006| |In- 

oue Fukui||2013| |Dobbs, Pringle Duarte-Cabral||2015{ ~ 

Moreover, instances where on-going star formation may have 
been triggered by cloud-cloud collisions are hard to identify, 
and the observations are difficult to interpret unambiguously 


(Haworth et al. 2015). Possible instances include Wester- 
lund 2 ( Furukawa et al. 1 2009 Ohama et al.|2010 ), the Trffid 


Nebula (Torii et al. 2011), and the Serpens Main Cluster 


( Duarte-Cabral et ai.||20li ). 


In this paper we present simulations of cloud-cloud col¬ 
lisions, and discuss the statistical properties and spatial dis¬ 
tributions of the protostars that form. We limit considera¬ 
tion to head-on collisions involving a single cloud mass and 
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a single cloud radius, and focus our attention on the con¬ 
sequences of varying the velocity of the collision, v^. The 
consequences of varying the cloud mass, the cloud radius, 
and the impact parameter of the collision, will be explored 
in subsequent papers. 

In Section!^ we describe the simulation set-up. In Sec¬ 
tion we present and discuss the results. The main conclu¬ 
sions are summarised in Section [4] 


2 SIMULATION SET-UP 
2.1 Initial conditions 

Initially, the two colliding clouds both have mass Mq = 
500 Mq, radius = 2 pc, uniform density = 1.01 x 
10“^^gcm“^, and hence freefall time tpp = 2.1 Myr. The 
fraction of hydrogen, by mass, is X = 0.70, and the hy¬ 
drogen is assumed to be molecular, so the initial molecular 
hydrogen number-density is = 2.1 x 10^H2cm“^. The 
gas is initially isothermal at Tq = 10 K, corresponding to an 
isothermal sound speed = 0.19 km s“^. The cloud centres 
are initially at 

{x,y,z) = 


(±2.02,0,0) pc. 


( 2 . 1 ) 
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i.e. the clouds are almost touching, and the initial bulk- 
velocities of the clouds are 


iyx-) Vy, Vz) — (T'^o 5 0? 0) • 


( 2 . 2 ) 


We presume that the two colliding clouds are part of 
a larger, approximately virialised cloud complex, with 
mass between 4 Mq = 2,000Mq and 40 Mq = 20, 000 


- an d hence, using Larson’s scaling relation (Larson 
1981), velocity dispersion between ~ 1.2kms“^ and ~ 
2.0kms“^. Consequently the clouds are given velocities 
Vq = 1.2, 1.4, 1.6, 1.8 and 2.0kms“^, and they collide 
with relative velocities in the range Auq ~ 2.4kms~^ to 


Av^ ~4.0 km s . We discuss at the end of Section 


3.1 


what 


happens if Av^ falls outside this range. The Mach Numbers 
of the accretion shocks bounding the shock-compressed layer 
are of order 6 to 10, giving compression factors in the range 
36 to 100. 


2.2 Constitutive physics 

The gas is evolved with a barotropic equation of state. 


(2.3) 


with Tq = lOK and Pcrit — 10“^^ gcm“^. At low-densities, 
P ^ PcRiT? is approximately isothermal at ~ 10 K. 

At high densities, p ^ Pcrit ? g^® approximately adi¬ 

abatic with adiabatic index 7^5/3 (corresponding to cool 
molecular hydrogen in which the rotational degrees of free¬ 
dom are not significantly excited, plus some monatomic he¬ 
lium) . This prescription mimics approximately the variation 
of temperature with density at the centre of a 1 Mq proto¬ 
star ( Masunaga Inutsuka||1999 ). 


2.3 Numerical method 

The simulations are performed with the Smoothed Particle 
Hydrodynamics code GANDALF (Rubber et ah, in prepara¬ 
tion). A binary tree is used to compute gravitational forces 
and collate particle neighbour lists. The grad-/i formulation 
of the evolution equations is invoked [Price Monaghan] 
(2004), and they are solved with a second-order leapfrog 


integrator. Artificial viscosity is treated using the method 


of Morris & Monaghan (1997), with parameters ol^.^ — 0.1, 
a — 1 and /3 = 2. Sinks are introduced and evolved using 
the standard procedure described in [Rubber, Walch Whit-] 
worth (2013), with Psink — 10“^^gcm“^; hence condensa- 
tions which are converted into sinks are already well into 
their Kelvin-Helmholtz contraction phase, and their proper¬ 
ties are essentially independent of the choice of Psink • 

All simulations are performed with A/"rpQrp ~ 10® SPH 
particles, so each SPH particle has mass ce 0.001 Mq . 

The resolution-parameter is set to p = 1.2, so that an SPH 
particle typically has ~57 neighbours. Hence the mass res¬ 
olution is ~ 0.1 Mq (i.e. just above the hydrogen-burning 
limit). The state of the simulation is output every 10, 000 yr. 

Feedback from protostars is not included, and therefore 
the simulations are terminated once 10% of the total mass 
has been converted into protostars, on the assumption that 
feedback would be terminating star formation around this 
stage. To include feedback would require us (a) to reduce 


the resolution (increased with the consequence that 

low-mass H-burning stars could not form, and (b) to intro¬ 
duce sinks at a lower density threshold (lower Psink)? ^o that 
the simulations could be followed further and massive stars 
(> 8 Mq) could form (cf. Dale, Ercolano & Bonnell[[2015 ). 
We are presently repeating the simulations presented here 
with larger and lower Pgi^K? with a view to under¬ 

standing the formation and evolution of bipolar Hii regions 
(e.g. Deharveng et al.[[2015 ), and these simulations will be 
presented in a future paper. 


2.4 Perturbations 

For each velocity, , we perform ten simulations. Of these 
ten simulations, five are performed with clouds that ini¬ 
tially have zero internal velocity dispersion (Set 1), and five 
are performed with clouds that initially have a very low- 
amplitude internal turbulent velocity field (Set 2). 

For the Set 1 simulations, the SPH particles are given 
random positions in a cube with side = 4 pc. The 

cube is then settled, using periodic boundary conditions, 
until the fractional density fluctuations have standard de¬ 
viation less than 0.4%. Finally, particles are culled if they 
fall outside a sphere of radius 2 pc. This sphere, and a clone 
of it rotated through random Euler angles, then represent 
the two clouds. The particles in a cloud have zero velocity 
dispersion. Initial density perturbations derive solely from 
the very small residual fluctuations in the particle density. 
The ratio of thermal to gravitational energy in an individual 
cloud is c^thermal =0-083 < 0.5, so the pre-collision clouds 
are not virialised. 

For the Set 2 simulations, the particles are distributed 
in the same way as Set 1. Next a random Gaussian tur¬ 
bulent velocity field, with power-spectrum Pk oc is 

generated on a 128® grid, and the velocities of the SPH 
particles are computed by interpolation on this grid. Fi¬ 
nally, the particle velocities are shifted to the centre of mass 
frame, and re-scaled so that they have a one-dimensional 
velocity dispersion ct^urb ~ 0.066 kms“^, corresponding 
to Mach Number Af^uRB — 0.35. The ratio of turbu¬ 
lent to gravitational energy in an individual cloud is then 
q^turb =0.010, so the pre-collision clouds are not virialised. 


THERM 4“ ^TURB 0.093 < 0.5. 


The initial conditions are not intended to be realistic. 
Realistic clouds would not have uniform density, and they 
would have a much higher level of turbulence. The initial 
conditions are chosen because we want to isolate the effects 
of the collision from the effects of pre-existing density or 
velocity sub-structures. Also, we want to keep the parameter 
space as small as possible, in order to be able to quantify 
cause and effect. We will explore in a future publication the 
extent to which the coherent large-scale structures generated 
here survive when the individual clouds have pre-existing 
density and velocity sub-structures of different amplitudes. 

If evolved in isolation, the individual clouds collapse and 
fragment to form protostars, but they do so on a somewhat 
longer timescale than we consider here ( ^ 2 Myr). 
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Figure 1. False-colour images of colliding clouds, as seen looking across the collision axis, at (i.e. when the first protostar forms) 
and (i.e. when 10% of the mass has been converted into protostars), (a) = 1.2kms“^, at = 1.33 Myr; (b) Vq = 1.2kms“^, at 

^ 10 % ~ 1.47Myr; (c) Vq = 2.0kms“^, at = 1.14Myr; (d) Vq = 2.0kms“^, at = 1.28Myr. Colour represents log^^ {S/gcm“^}, 
where S is column-density projected on the (x, y)-plane. The black dots are protostars. 


3 RESULTS 
3.1 Overview 

The cloud-cloud collision creates a dense shock-compressed 
layer on either side of the x = 0 plane, as illustrated in Fig. 

Once this layer becomes sufficiently massive, it fragments 
gravitationally to produce a network of filaments that feed 
material into dense prestellar cores, and hence into proto¬ 
stars. At the same time the layer starts to collapse laterally, 
i.e. towards the x-axis. 

The statistical properties of the protostars that form, 
and their initial spatial distribution, do not depend on how 
the initial perturbations are seeded. In other words there 
is no discernible difference between the simulations of Set 
1 (small density perturbations) and those of Set 2 (small 
density and velocity perturbations). Therefore we merge the 
results from the two sets, and discuss them as a single set. 

The statistical properties of the protostars that form, 
and their initial spatial distribution, do depend, critically, 
on the velocity of the collision, Uq, and we therefore focus 
on this dependence. We identify two critical times, is 
the time at which the first protostar forms (i.e. the first sink 
particle); Fig. presents column-density images at for 
collisions involving the same initial clouds (same initial par¬ 
ticle positions and same initial turbulent velocity field) but 
different bulk velocities, v^. is the time at which the 


total mass of protostars first exceeds 100 Mq , in other words 
10% of the mass has been converted into protostars; Fig. 
presents column-density images at , for the same simu¬ 
lations presented in Fig. All simulations are terminated 
at on the assumption that the inclusion of feedback 

would decelerate, or even terminate, star formation around 
this stage - and also because, by this stage, the progress 
of a simulation has normally slowed so much that it is not 
practical to compute much further. 

For lower , the surface-density of the layer builds up 
more slowly, and therefore there is more time for lateral 
contraction of the layer (i.e. shrinkage towards the colli¬ 
sion axis), before it becomes gravitationally unstable against 
fragmentation. This leads to a network of filaments that are 
predominantly radial, and therefore appear like spokes, feed¬ 
ing material into a single, dominant, central hub. Some of 
the material arriving in the hub has already condensed out 
as protostars, but some remains as interstellar gas. Conse¬ 
quently, the protostars in the hub compete for this inflowing 
gas, and one or two grow to very high mass. At the same 
time, many are ejected by dynamical interactions. By 
there is an extensive diaspora of low-mass protostars, around 
a single large star cluster. 

For higher Uq, the surface-density of the layer builds 
up more quickly, and so it becomes gravitationally unsta¬ 
ble against fragmentation more quickly. As a consequence. 
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Figure 2. False-colour images of fragmenting layers at the moment the first protostar forms, , from simulations using the same initial 
clouds (one of the ones that start with weak density and velocity perturbations) but different bulk velocities, Vq. (a) =1.2kms“^, 

= 1.33Myr. (b) = 1.4kms“^, = 1.28Myr. (c) = 1.6kms“^, = 1.23Myr. (d) = 1.8kms“^, = 1.19Myr. (e) 

Uq = 2.0kms“^, = 1.14Myr. Colour represents log^^ {E/gcm“^}, where S is column-density projected on the (y, 2 ;)-plane, i.e. 

looking along the direction of the collision. The first protostar is represented by a black dot. 


the layer has less time to contract laterally, and therefore 
it fragments into a network of filaments that is more like a 
spider’s web, with many intersections. Prestellar cores con¬ 
dense out at each of these intersections, but, since these 
cores can only accrete a limited amount of material (i.e. 
just from the the small sections of filament closest to them), 
each core spawns only a small sub-cluster, typically compris¬ 
ing one moderately massive protostar and a few lower-mass 
protostars. Consequently, by there have been fewer dy¬ 
namical ejections, and the individual sub-clusters are still 
relatively isolated from one another. We anticipate that if 
the evolution were followed further, these sub-clusters would 
fall together and merge. 


With lower v^, the maximum stellar mass is larger 
(than with higher but the mean mass, and the minimum 
mass are both lower. With lower , the instantaneous pro- 
tostellar accretion rates show a greater spread (than with 
higher Vq), even at high masses, reflecting the stochastic 
nature of competitive accretion. 



Figure 4. The time of formation of the first protostar, (red 
stars, red best-fit line), and the time at which 10% of the total 
mass has been converted into protostars, (blue filled circles, 
blue best-fit line), plotted against the cloud velocity, Uq. The 
green dashed line is the value of predicted by Eqn. (|3.3|>. 


3.2 Fragmentation of the layer 

Fig. E shows the filamentary structure that has developed 
by the time the first protostar forms at , for simulations 


log column density g cm' 
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Figure 3. False-colour images of fragmenting layers at the moment, , when 10% of the mass has been converted into protostars, from 
the same simulations presented in Fig. I (a)'f;Q=1.2kms ^= 1.47 Myr. (b) Uq = 1.4 km s ^= 1.42 Myr. (c) Uq = 1.6 km s 
^ 10 % =l-37Myr. (d) =1.8kms“^, =1.33Myr. (e) =2.0kms“^, =1.28Myr. Colour represents log^^ {S/gcm“^}, where E 

is column-density projected on the (y, zj-plane, i.e. looking along the direction of the collision. The black dots represent the protostars. 


with different collision velocities, Vq. In all cases the same 
initial cloud has been used, in this case a cloud with uniform 
density and a weak turbulent velocity field. The turbulent 
velocity field, and the cloud orientation, are the same in all 
five simulations, and consequently the configuration of the 
filaments at is broadly similar. Indeed, there is a clear 
predominance of filaments with loci y ~ — 1.1 z, reflect¬ 

ing the fact that in this case the initial turbulent velocity 
field has a significant fraction of its power in a compres- 
sional mode with ~ +l.l/Cy. This indicates that filaments 
growth is seeded largely by the initial perturbations, even 
though these are small (see Section 2.4). 

The only difference is that with lower collision velocity 
the surface-density of the layer builds up more slowly, and 
so the filamentary structure takes longer to develop; as a 
result, by the time the first protostar forms, lateral contrac¬ 
tion of the layer has progressed further, and the extent of 
the filamentary network is smaller. 


3.3 Formation of the first protostar 

A shock-compressed layer created by two semi-infinite anti¬ 
parallel flows with density and velocity fragments at 


time 


“ (23/2Gpo?;o) 

(Whitworth 2015, in preparation). Here is the effec¬ 
tive isothermal sound speed in the shocked gas, and ‘semi¬ 
infinite’ means that, if the contact discontinuity between the 
flows is the x = 0 plane, the gas has infinite extent in the y- 
and ^-directions. 

However, the layer created by a cloud-cloud collision has 
finite extent in the y- and ; 2 ;-directions, and consequently it 
contracts towards the x-axis on a time-scale 


(32GP0) ’ 


(3.2) 


thereby accelerating the condensation of the first protostar. 

We therefore posit - somewhat arbitrarily - that the 
time at which the first protostar forms might be approxi¬ 
mated by the logarithmic mean of and tpp, i.e. 


— (^crit^ff) 

1 


1/2 


3 TT ttg 
2 ^ 7 ^ 


1/4 


(3.3) 


2(Gpo)^/^ 


log column density / g cm. 
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Setting ttg = 0.19kms 


^ this reduces to 

- 1/4 


(3.4) 


Fig.|^shows the values of and from the simulations, 
plotted against iIq, and the predictions of Eqn. (3.3). 


3.4 Reconfiguration of the filaments 

After , the network of filaments continues to be amplified, 
and - in the cases illustrated in Figs. an di the predomi¬ 
nance of filaments with ce -f-l.lky wanes, as other modes 
in the initial turbulent velocity field assert themselves, to 
produce a criss-crossing network of filaments. At the same 
time, the shock-compressed layer contracts laterally. This 
contraction is non-homologous, because the surface-density 
of the layer increases towards the x-axis, and consequently 
the inner parts of the layer converge on the centre on a 
shorter timescale than the outer parts. This has the effect of 
stretching filaments into predominantly radial orientations. 

At lower collision velocities, Xq < 1.6 kms“^, radial 
stretching is the main process reconhguring the filaments. 
This has two consequences. First, the filaments morph into 
a system of radial arteries feeding matter into the centre, i.e. 
a hub and spokes system (see the frames on the top row of 
Fig.§. Second, wherever a filament fragments to produce 
protostars, these protostars are diverging from one another, 
due to the tidal stretching of the filament out of which they 
have just condensed. Therefore they tend to have low masses 
and they tend not to interact with one another until they 
reach the central hub. There the protostars interact violently 
with one another, and most of the low-mass ones are ejected, 
to produce a diaspora of low-mass protostars, whilst a few 
settle into the centre and grow by competitive accretion. 

At higher collision velocities, > 1.6kms“^, filaments 
grow from many different turbulent modes, producing a 
criss-crossing network, with multiple nodes wherever two or 
more filaments intersect (see the frames on the bottom row 
of Fig. 1^. Each node develops into a prestellar core, but 
since the amount of matter in the short hlamentary seg¬ 
ments between nodes is small and is shared between two 
nodes, these cores normally have quite low mass. Conse¬ 
quently, each individual core can only spawn a small sub¬ 
cluster of protostars, and typically this sub-cluster contains 
just one moderately massive protostar and a few less massive 
ones. 



1.6 1.8 
Velocity of clouds / n 


Figure 5. Protostellar mass ranges. For each collision velocity, 
Vq , we combine the data from all ten realisations. The filled cir¬ 
cles give the median mass. Scanning downwards (upwards) from 
the median mass, the first blue line gives the 25^^ (75^^) centile, 
i.e. the interquartile range; the next red line gives the 10^^ (90^^) 
centile; and the last green line gives the 0^^ (100^^) centile (i.e. 
the full mass range). 



Figure 6. The cumulative mass of protostars at ^ func¬ 

tion of protostellar mass: red curve, Vq =1.2kms“^; blue curve, 
Vq = 1.4kms“^; green curve, = 1.6kms“^; black curve, 
Vq =1.8kms“^; purple curve, Vq =2.0kms“^. 


At lower collision velocity, star formation starts later, but 
there is less mass in the shock-compressed layer, when it 
fragments. This is compensated by the fact that the layer has 
had more time to contract laterally, so the surface-density of 
the layer is comparable, and hence the star formation rate 
is similar. 


3.5 The rate of star formation 


The simulations are terminated at when 10% of the 

total mass, i.e. 100 Mq, has been converted into protostars. 
Values of are plotted against on Fig.j^ The duration 
of star formation. 


Atgp = ^ 10 % ~ ^*1 — 0.14 ± 0.01 Myr, (3-5) 


is essentially independent of the collision velocity, v^. Thus, 
the overall star formation rate in the layer is also essentially 
independent of , 


100 Mq 


(3.6) 


Table 1. Parameters of the mass functions of protostars formed 
with different collision velocities. Column 1 gives the collision ve¬ 
locity, Uq, column 2 gives the mean number of protostars, A/"^, 
formed from one collision (and totalling 100 Mq), and columns 
3 through 6 give, respectively, the mean, , standard devi¬ 
ation, , skew, , and kurtosis, , of the protostellar 

log-masses, /Mq) 


'^o 

K 

l^eoqM 

^(oqM 

'^(oqM 

'^(oqM 

1.2 kms“^ 

198 

-0.65 

+0.45 

+1.92 

+4.08 

1.4kms“^ 

163 

-0.56 

+0.53 

+1.41 

+1.12 

1.6 kms“^ 

108 

-0.29 

+0.60 

+0.43 

-1.21 

1.8kms-i 

88 

-0.13 

+0.58 

-0.11 

-1.33 

2.0kms“^ 

86 

-0.01 

+0.50 

-0.58 

-0.72 


M, 


700 Mq Myr“^ . 
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3.6 The mass function of the protostars 

As the collision velocity, is increased, the median proto- 
stellar mass increases, but the maximum protostellar mass 
and the spread of protostellar masses both decrease. Lower 
Vq leads to the formation of a single central hub, into which 
the majority of protostars fall; there they partake in com¬ 
petitive accretion, and consequently a small number of pro¬ 
tostars grow to very large masses, and a large number of 
protostars are ejected dynamically, producing an extensive 
diaspora of predominantly low-mass protostars. At higher 
Vq, there are many different cores spawning protostars in¬ 
dependently of one another. In each such core a moderately 
massive protostar forms, along with a few somewhat lower- 
mass protostars; in these small-iV sub-clusters, dynamical 
ejections do not occur so quickly, and when the simulations 
are terminated at most protostars are still near the lo¬ 
cation where they were born. 

Fig. E illustrates these trends by plotting the mass 
ranges, as a function of collision velocity. The green lines 
represent the complete range of masses at each velocity. The 
upper green line here is significant, and reflects the fact that 
the most massive protostars are formed in the lower-r^Q col¬ 
lisions. The lower green line simply reflects a single very 
low-mass protostar formed in one of the ten simulations at 
that particular collision velocity, and is therefore less signifi¬ 
cant. The red lines represent the range between the 10^^ and 
90^^ centiles; in other words, by number^ 10% of stars fall be¬ 
low this line, and 10% fall above it. Similarly, the blue lines 
represent the range between the 25^^ and 75^^ centiles, also 
known as the interquartile range. The filled circles represent 
the median masses for each collision velocity. 

Fig. shows the cumulative mass of protostars at 
as a function of protostellar mass, at the different collision 
velocities. This illustrates how the number of low-mass pro¬ 
tostars increases dramatically with decreasing collision ve¬ 
locity, i’q, once Vq falls below ~ 1.6kms“^. It also illus¬ 
trates how the two most massive protostars become more 
massive with decreasing as falls below ~ 1.6kms“^. 
For higher collision velocities, > 1.6kms“^, the distribu¬ 
tion of masses is much tighter, i.e. it has a smaller standard 
deviation (see Table [^. 

In no case is the distribution of masses log-normal. Ta¬ 
ble gives the mean, standard deviation, cr^^, skew, 

7 &gM, and kurtosis, of log^jj(M*/M 0 ), for different Vq. 

As Vq increases, increases monotonically, i.e. the mean 
stellar mass increases; (7^^ increases and then decreases, 
but does not vary by much (0.45 < < 0.60); 7^^ de¬ 

creases monotonically, indicating that the mass function is 
very assymetric at low (due to the large number of very 
low-mass protostars), but becomes increasingly symmetric 
at larger ; the variation of the kurtosis indicates that the 
mass function at low is much more sharply peaked than 
at high Vq . 

Although on large scales the observed Initial Mass Func¬ 
tion (IMF) appears to vary very little, it is still possible 
that there exists some scale below which there are signifi¬ 
cant variations, over and above those that arise purely from 
sampling statistics; the observed universal IMF then con¬ 
stitutes an average over a representative ensemble of such 
regions. The results above are a prediction for what might 
be uncovered if the stellar population resulting from a sin¬ 


gle cloud-cloud collision could be accurately collated, at an 
early stage, before it is dispersed and/or contaminated with 
other stars. 


3.7 Accretion rates 

Fig.j^shows the final protostellar masses, M^, at plot¬ 
ted against their time of formation, tpoRM? ^^e time at 
which the corresponding sink particle was created. As the 
collision velocity is increased, the formation time time de¬ 
creases, the maximum mass decreases, and the number of 
low-mass protostars decreases. 

Fig. shows the protostellar accretion rates, M^, plot¬ 
ted against protostellar mass, M^, both evaluated at the 
end of the simulation, For lower collision velocities, 

Vq < 1.6kms“^, final masses are concentrated towards low 
values, 0.08 to 3 Mq , and there is a wide range of final accre¬ 
tion rates at all masses, 10“^ to 10^ M^ Myr“^, apart from 
the most massive protostars, which are all still growing very 
rapidly. For higher collision velocities, Vq > 1.6kms“^, final 
masses and final accretion rates are concentrated at quite 
high values, 1 to 10 M^ and 10 to 100 M^ Myr“^, respec¬ 
tively. 

3.8 Spatial distribution of protostars 

We have analysed the spatial distributions of the protostars 
formed in the simulations, and how they evolve, by project¬ 
ing the initial and final positions of protostars onto the x^O 
plane and constructing Minimal Spanning Trees (MSTs). 
The top row of Fig.j^shows an example of the MSTs for the 
initial (left) and hnal (right) positions obtained in a low- 
velocity collision, with = 1.2kms“^. The bottom row of 
Fig. IE shows an example of the MSTs for the initial (left) 
and final (right) positions obtained in a high-velocity colli¬ 
sion, with Vq = 2.0kms“^. We see that in the aftermath of 
a low-velocity collision the initial positions of the protostars 
(Fig. [^) are strongly concentrated in a central cluster, but 
the final positions (Fig. Ep) are spread out over a large area, 
due to violent dynamical ejections of low-mass protostars. 
In contrast, in the aftermath of a high-velocity collision the 
initial positions (Fig.|^) are quite spread out, and the final 
positions (Fig. Et) are, statistically, not much changed, re¬ 
flecting the fact that thus far there has been little dynamical 
evolution within the separate sub-clusters, and the proto¬ 
stars are still quite close to the places where they formed. 

This is quantified in Fig. where we plot the mean 
edge-lengths on the initial and final MSTs against the colli¬ 
sion velocity, . The mean edge-length on the initial MST 
is a strong function of the collision velocity, increasing by 
almost an order of magnitude (lower blue curve), as we 
go from Vq = 1.2kms“^ (very concentrated single cluster) 
to Vq — 2.0kms“^ (many more-or-less independent sub¬ 
clusters). The mean edge-length on the final MST is almost 
independent of the collision velocity (upper green curve). 
For lower-velocity collisions, this is largely due to the dis¬ 
persal of the low-mass protostars ejected from the central 
cluster. As the collision velocity is increased, this effect be¬ 
comes increasingly unimportant, and the mean edge-lengths 
on the initial and final MSTs obtained with =2.0kms“^ 
are indistinguishable. 
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Figure 7. The final mass of each protostar, at (the end of the simulation) plotted against its time of formation, tpoRM^ 

the time at which the corresponding sink particle is created. Red dots, = 1.2kms“^; blue crosses, Vq = 1.4kms“^; green stars, 
Vq =1.6kms“^; black open circles, Vq =1.8kms“^; purple triangles, =2.0kms“^. 


In order to estimate the Q-parameter, and hence quan¬ 
tify further the evolution of the clusters, we must normalise 
the edge-lengths on the MST according to 

^ {"^mstI > (3-7) 

^ ' MST 

where is the total number of protostars and is 

one of the edges on the MST (see Cartwright & Whitworth 
2004). In addition, we must evaluate the mean of the edges 
connecting all pairs of protostars, 

z=AA^-1 j=J\f^ 

’ - ww^) ^ ^ 

* ^ ' 1=1 j=t+l 


where Sij is the edge between protostars i and j. Then the 
Q-parameter is given by 


Q = ^. 

S 


(3.9) 


Fig. 11 shows the mean location on the (m, s)-plane 


for the clusters produced with different v^. For each value 
of i’q, the symbol representing the initial positions of the 
protostars is smaller than the symbol representing the fi¬ 
nal positions of the protostars. The dashed lines are lines of 
constant Q. In all cases, the cluster evolves towards smaller 
rfi and smaller s, but larger Q. This reflects two factors. 
First, locally protostars are becoming more tightly bound 
to their nearest neighbours, which makes the mean edge- 
lengths (both ffi and s) smaller. Second, the clusters formed 
following lower-velocity collisions are strongly centrally con¬ 
densed, and become even more so, as the higher-mass proto¬ 
stars settle into the centre, and the lower-mass protostars are 


ejected at high speed; conversely, the clusters formed follow¬ 
ing higher-velocity collisions are made up of many more-or- 
less independent sub-clusters, formed at the intersections of 
a spider’s web of filaments, and at the end these sub-clusters 
are only just beginning to dissolve. 


3.9 More extreme collision velocities 

If the clouds considered here (500 Mq , 2pc) have lower veloc¬ 
ities than we have discussed, i’q < 1 kms“^, they collapse in¬ 
dividually (and form separate star clusters) faster than they 
collide, so the effect of the collision is simply to merge two ex¬ 
isting star clusters. Conversely, if the clouds have higher ve¬ 
locities than we have discussed, >2.4kms“^, the trends 
described above continue: a more extended spider’s web of 
intersecting filaments forms; the spacing of the filaments de¬ 
creases, and hence the masses of the subclusters forming at 
their intersections decreases; star formation occurs earlier; 
the median mass tends to a constant value; and the range 
of masses continues to decrease. We have tested this trend 
up to Vq — 4kms“^, At'q = 8 kms“^. 


4 CONCLUSIONS 

We have explored the phenomenology of star formation 
triggered by cloud/cloud collisions, specifically what hap¬ 
pens when two quiescent and approximately uniform-density 
clouds with mass 500 Mq and radius 2 pc collide head-on at 
relative speeds between 2.4 and 4.0kms“^. 

(i) In all cases the colliding clouds produce a shock- 
compressed layer, which fragments into filaments. 
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Figure 8. Log-log plot of the rate of accretion, against the 
protostellar mass, , both evaluated at the end of the simu¬ 
lation, Red dots, Vq = 1.2kms“^; blue crosses, = 

1.4kms“^; green stars, = 1.6kms“^; black open circles, 
= 1.8kms“^; purple triangles, = 2.0kms“^. Note that 
the lowest accretion rates are quantised at values 0.1 Mq Myr“^, 
0.2MQMyr“^, O-SM^Myr”^, etc., corresponding to sinks ac¬ 
creting one, two, three, etc. SPH particles (hence 10“^ Mq , 
2 X 10“^ Mq 3 X 10“^ Mq, etc.), during the 0.01 Myr interval 
between the last two outputs. Similarly some accretion rates fall 
on the line = M^/O.OlMyr, corresponding to sinks formed 
during the 0.01 Myr interval between the last two outputs. 





Figure 9. Minimal Spanning Trees connecting the positions of 
protostars projected on the x = 0 plane (the collision axis). The 
top row represents a simulation with = 1.2kms“^; (a) shows 
the initial positions (i.e. where the protostars formed), and (b) 
shows the final positions (i.e. at tio%)- bottom row represents 
a simulation with Vq = 2.0kms“^; (c) shows the initial positions, 
and (d) the final positions. 



Collision velocity 


Figure 10. The mean edge-lengths from the initial and final 
MSTs, as a function of collision velocity, . The lower blue curve 
represents the initial MSTs, and the upper green curve represents 
the final MSTs. 


(ii) At the lower relative velocities, it takes longer before 
fragmentation becomes non-linear and delivers protostars. 

(iii) At the lower relative velocities the lateral collapse 
of the layer drags and stretches the filaments towards the 
collision axis, so that they become predominantly radial. 
Protostars that condense out of the filaments fall along the 
filaments into the centre, along with residual gas, and a sin¬ 
gle massive cluster forms there. Competitive accretion leads 
to the growth of a few very massive protostars, and the 
ejection of a large number of low-mass protostars. The dis¬ 
tribution of filaments around the single massive cluster can 
be described as a hub and spokes system. 

(iv) At the higher relative velocities, there is less time 
for lateral collapse of the layer, and so the filaments form 
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Figure 11. The locus of simulations on the (m, s)-plane, where m 
and s are defined in Eqns. and ( |3.8| > . Red dots mark the mean 
for Uq = 1.2kms“^; blue crosses, for Vq = 1.4kms“^; green 
stars, for Uq = 1.6 km s“^; black open circles, for Uq =1.8 kms“^; 
and purple triangles, for Uq =2.0kms“^. For each Uq the symbol 
is smaller for the initial positions and large for the final positions. 
The dashed lines are lines of constant Q = fh/s (Cartwright & 
Whitworth 2004). 


a network like a spider’s web, with many more-or-less inde¬ 
pendent star-forming cores collecting at the intersections of 
the filaments. Since these cores grow from relatively short 
filamentary segments, they have limited mass, and typically 
spawn only a small number of protostars, usually one fairly 
massive one and a few smaller ones. Because there is less 
mass, there are fewer protostars in each core, so dynamical 
ejections occur later and are less violent. 

(v) As the collision velocity is increased, the mean proto- 
stellar mass increases, but the range of protostellar masses 
goes down, so the maximum mass actually decreases and 
the minimum mass increases; there are many fewer low-mass 
protostars. 

(vi) The spatial distribution of protostars when they first 
form is marginally fractal, due to the underlying filamen¬ 
tary distribution of the gas. At lower collision velocities, the 
protostars quickly relax to a single very centrally condensed 
configuration. At higher collision velocities, the protostars 
are still in small independent sub-clusters, close to where 
they formed, at the end of the simulation. 
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